PROCEEDINGS 

OF SCIENCE 



o 

C 



in: 



> 

in 
o 



X 



Critical properties of the two-dimensional Z(5) 
vector model 



Oleg Borisenko 

Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine 
03680 Kiev, Ukraine 



E-mail: oleg@bitp.kiev.ua 



Gennaro Cortese* 

Dipartimento di Fisica, Universita delta Calabria, and INFN - Gruppo Collegato di Cosenza 
1-87036 Rende, Italy 



E-mail: cortese@cs . inf n . it 



Roberto Fiore 

Dipartimento di Fisica, Universita delta Calabria, and INFN - Gruppo Collegato di Cosenza 
1-87036 Rende, Italy 
E-mail: 



fiore@cs.infn.it 



Mario Gravina 

Laboratoire de Physique Theorique, Universite de Paris-Sud 11, Batiment 210 

91405 Orsay Cedex, France 

E-mail: |Mario . Gravina@th . u-psud. f r 



Alessandro Papa 

Dipartimento di Fisica, Universita della Calabria, and INFN - Gruppo Collegato di Cosenza 
1-87036 Rende, Italy 



E-mail: papa@cs . inf n . it 



The two-dimensional Z(5) vector model is investigated through the determination of critical points 
and one critical index. To this purpose a new cluster algorithm has been developed valid for 2D 
Z(N) models with odd values of N. Results are compared with analytical predictions. 



The XXVIII International Symposium on Lattice Field Theory, Lattice2010 
June 14-19, 2010 
Villasimius, Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Critical properties of the 2D Z(5) vector model 



Gennaro Cortese 



1. Introduction 



The Berezinskii-Kosterlitz-Thouless (BKT) phase transition is known to take place in a variety 
of two-dimensional (2D) systems, the most common being the 2D XY model [j]}. Here we are going 
to study an example of lattice spin model where this type of the transition shows up, namely the 2D 
Z(N) spin model, also known as vector Potts model. On a 2D lattice A = L 2 with linear extension 
L and periodic boundary conditions, the partition function of the model can be written as 
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(1.1) 

in the standard formulation with N — 1 different couplings. 

Some details of the critical behavior of 2D Z(N) spin models are well known - see the review in 
Ref. The Z(/V) spin model in the Villain formulation has been studied analytically in Refs. [g]. 
It was shown that the model has at least two phase transitions when N > 5. The intermediate phase 
is a massless phase with power-like decay of the correlation function. The critical index tj has been 
estimated both from the renormalization group (RG) approach of the Kosterlitz-Thouless type and 
from the weak-coupling series for susceptibility. It turns out that T](/3i^) = 1/4 at the transition 
point from the strong coupling (high-temperature) phase to the massless phase, i.e. the behavior is 
similar to that of the XY model. At the transition point jS^ 2 ' from the massless phase to the ordered 
low-temperature phase one has 17 (j3 c ) = 4/N 2 . A rigorous proof that the BKT phase transition 
does take place, and so that the massless phase exists, has been constructed in Ref. [Q] for both 
Villain and standard formulations (with one non- vanishing coupling j8i). Monte-Carlo simulations 
of the standard version with A 7 = 6, 8, 12 were performed in Ref. [^]. Results for the critical index 
tj agree well with the analytical predictions obtained from the Villain formulation of the model. 

Here we investigate the case N = 5, the lowest number where the BKT transition is expected. 
The motivation of our study is two-fold: (i) to compute critical indices at the transition points, 
which could serve as checking point of universality; (ii) to develop and test a Monte Carlo cluster 
algorithm valid for odd values of N, not yet available in the literature, to our knowledge. 



2. Algorithm and numerical set-up 

In this work we concentrate our attention to the model defined by Eq. ( |1.1[ ) with only one 
non-zero coupling, j8i = /3 . The Hamiltonian of the model is 

H = -p^cos (^(si - sj)) , Si = 0,l,...,N-l, (2.1) 

(ij) V iy J 

with summation taken over nearest-neighbor sites. For N = 2 this is the Ising model, whereas in 
the N — > 00 limit we get the XY model. 

A cluster algorithm for the Monte Carlo numerical simulation of this model is available in the 
literature only for even N [Bt]. Here we develop a new algorithm, valid instead for odd Af, by which 
an accurate numerical study of the model can be performed for N = 5, i.e. the smallest Af value 
for which the phase structure described in the Introduction holds. Here are the steps of our cluster 
algorithm for the update of a spin configuration {st}: 
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Figure 1: Scatter plot of the complex magnetization Ml at j3 — 0.80, 1.10, 1.50 in Z(5) on a 64 2 lattice. 

• choose randomly n in the set {0, 1,2, ... ,N — 1} 

• build a cluster configuration according to the following probability of bond activation be- 
tween neighboring sites ij 



I 1 -exp(-26 a,a,) if a,a, > ., . (2%, 

'H otherwise ' wxth = sm 



• "flip" each cluster, with probability 1/2, by replacing all spins belonging to it according to 
the transformation 

Si ->■ mod(-5; + 2n + N,N) , 
which amounts to replacing each spin si in a cluster by the spin sj for which (Xj = —a,. 

It is easy to prove that this cluster algorithm fulfills the detailed balance. We have tested the effi- 
ciency of the cluster algorithm against the standard heat-bath algorithm and found that the cluster 
algorithm is strongly preferable (see Ref. [^] for details). 

The three phases exhibited by the 2D Z(5) spin model can be characterized by means of two 
observables: the complex magnetization Mi and the population Sl- The complex magnetization is 
given by 

m l = 4 L ex p = |Mz - |e ' v • (2 - 2) 



v- 

In Fig. [j] we show the scatter plot of Mi on a lattice with L = 64 in Z(5) at three values of j8, each 
representative of a different phase: j8 = 0.80 (high-temperature, disordered phase), j8 = 1.10 (BKT 
massless phase) and j8 = 1.50 (low-temperature, ordered phase). As we can see we pass from a 
uniform distribution (low j8) to a ring distribution (intermediate j3) and finally to five isolated spots 
(high j6). The naive average of the complex magnetization gives constantly zero, therefore Ml is 
not an order parameter. A convenient observable to detect the transition from one phase to the 
other is instead the absolute value \Ml \ of the complex magnetization. In Fig. || (left) we show the 
behavior of the susceptibility of \Ml \ , 



-L 2 ((\M L \ 2 ) 



(\M L \Y 



(2.3) 
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in Z (5) on lattices with L ranging from 16 to 1024 over a wide interval of j3 values. On each lattice 
x| M ' clearly exhibits two peaks, the first of them, more pronounced than the other, identifies the 
pseudo-critical coupling j3pd (L) at which the transition from the disordered to the massless phase 
occurs, whereas the second corresponds to the pseudo-critical coupling j8p C (L) of the transition 
from the massless to the ordered phase. It is evident from Fig. || that \Ml\ is particularly sensitive 
to the first transition, thus making this observable the best candidate for studying its properties. 

As a local order parameter to better detect the second transition, i.e. that from the massless to 
the ordered phase, we chose instead the population Sl, defined as 



N 



N-l 



max, =a jv_iQ;) 1 
L 2 N 



(2.4) 



where represents the number of spins of a given configuration which are in the state s/. In a phase 
in which there is not a preferred spin direction in the system (disorder), we have tij ~ L 2 /N for each 
index i, therefore Sl ~ 0. Otherwise, in a phase in which there is a preferred spin direction (order), 
we have rij ~ L 2 for a given index i, therefore Sl ~ 1 - In Fig. [|| (right) we show the behavior of the 
susceptibility of Sl, 

X { l S) =L 2 ((S 2 l )-(Sl) 2 ), (2.5) 

in Z(5) on lattices with L ranging from 16 to 1024 over a wide interval of j8 values. Again the peaks 
signalling the two transitions are clearly visible and their positions agree with Fig. ||, but now the 
second one is more pronounced. 

Other observables which have been used in this work are the following: 

• the real part of the "rotated" magnetization, Mr = \Ml\cos(5y) 

• the order parameter introduced in Ref. [|7|], niy = cos(5i//), 



where y/ is the phase of the complex magnetization defined in Eq. (2.2). For all observables con- 
sidered in this work we collected typically 100k measurements, on configurations separated by 
10 updating sweeps. For each new run the first 10k configurations were discarded to ensure ther- 
malization. Data analysis was performed by the jackknife method over bins at different blocking 
levels. 



3. Numerical results 

The first peak in the plot of the susceptibility x^ (see Fig. ^ (left)) indicates the transition 
from the disordered to the massless phase, while the second peak in the plot of the susceptibil- 
ity Xl^ ( see Fig- H^ght)) indicates the transition from the massless to the ordered phase. The 
couplings where these transitions occur (from now on denoted as the pseudo-critical couplings 

(l 2) 

j8p C ' (L)) have been determined by a Lorentzian interpolation around the peak of the correspond- 
ing susceptibility. Their values are summarized in Table [jj 

In order to apply the finite size scaling (FSS) program, the location of the infinite volume 
critical couplings j3 c ^ and /3j 2 ^ is needed. In Refs. [R ^] this was done by extrapolating the pseudo- 
critical couplings to the infinite volume limit, according to a suitable scaling law. First order 
transitions are ruled out by data in Table [jj. Second order transitions, though not incompatible 
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Figure 2: Behavior of the susceptibilities % L (left) and %l (right) versus j3 in Z(5) on lattices with several 
values of L. 

Table 1: Values of J3£P and in Z(5) on L 2 lattices. 



L 




As' 


16 


0.8523(20) 


1.1323(19) 


32 


0.91429(90) 


1.1363(11) 


64 


0.95373(40) 


1.13212(60) 


128 


0.98054(30) 


1.12875(66) 


256 


0.99838(20) 


1.12290(16) 


384 


1.00621(10) 


1.12103(50) 


512 


1.01112(20) 


1.11912(28) 


1024 


1.01991(10) 


1.11596(38) 



with data in Table |], are to be excluded, due to the vanishing of the long distance correlations 
combined with the clusterization property. Therefore, we assume that both transitions are of BKT 
type and adopt the scaling law dictated by the essential scaling of the BKT transition, i.e. 1~ ~ e ht v , 
which reads 

/3£' 2) =AP>+ Al ' 2 i ■ (3-D 
(\nL+B h2 )v 

The index v characterizes the universality class of the system. For example, v = 1/2 holds for the 
2D XY universality class. 

(I 2) 

Unfortunately, 4-parameter fits of the data for /3 pc ' (L) give very unstable results for the pa- 
rameters. This led us to move to 3-parameter fits of the data, with v fixed at 1/2. We found, as best 
fits with the MINUIT optimization code, 

/3 C (1) = 1.0602(20) A y = -2.09(20) B { = 0.27(18) £ 2 /d.o.f. = 0.48 L min = 64 
pP = 1.1042(12) A 2 = 0.578(41) B 2 = 0. Z 2 /d.o.f. = 0.61 L min = 128 

for the first and second transition, respectively. We observe that j8 c ^ is not far from the value 

(2) 

Of on the largest available lattice, thus supporting the reliability of the extrapolation to the 
thermodynamic limit. This is not the case for j3 c , suggesting that the considered volumes could 
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not be large enough for using the scaling law (CO). For this reason, we turned to an independent 

(12) 

method for the determination of j3 c , based on the use of Binder cumulants. 

In particular, for the study of the first transition, we considered the reduced 4-th order Binder 



cumulant l/[ M ^ defined as 



and the cumulant B^' defined as 



(M) (|M L | 4 ) 



{Mr) _ (\Mr — (Mr)\ 4 ) 

* 4 (|M«-(M«)|2\2' U " J 



while for the second transition we adopted again B^ R ^ and the cumulant B^ n ^ defined as 



Plots of the various Binder cumulants versus j8 show that data obtained on different lattice volumes 
align on curves that cross in two points, corresponding to the two transitions. We determined 
the crossing points by two methods: (i) by interpolating with polynomial lines data on different 
lattices near the crossing points and by looking for the intersection of these lines; (ii) by plotting 
the Binder cumulants versus (/3 — jS c )(logL) 1 / v , with v fixed at 1/2, and by looking for the optimal 
overlap of data from different lattices, by the % 2 method. As a result of this analysis (for details, 
see Ref. p) we arrived at the following estimates: /3 C (1) = 1.0510(10) and j8 c (2) = 1.1048(10). 
While pc is compatible with the infinite volume extrapolation of the corresponding pseudocritical 
couplings, /3 C is n °t, thus confirming the previous worries about the safety of the infinite volume 
extrapolation of j3pc • It should be noted, however, that a fit to j3pc (L) with the law ( |Q| ) and with 
the parameter j3 c fixed at 1.0510 gives a good # 2 /d.o.f., if only the three largest volumes are 
considered in the fit. 

The next step would be to extract other critical indices and check the hyperscaling relations 
at the two transitions. This calls for the FSS of magnetizations and susceptibilities at the critical 
couplings j3 c , which is in progress [Q]. We present here only two determinations of the effective 
v\ index, defined in Ref. [j|] as 

_ iog[r(j?)/r(j? )] 
^ {R) = iog[* /*] ' ( } 

where T(R) is the spin-spin correlation function and Ro an arbitrary parameter, chosen here equal 
to 10. This quantity is constructed in such a way that it exhibits a plateau in R if the correlator 
obeys the law 

m - , (3-6) 

valid in the BKT phase, j8 c (1) < j3 < j8 c (2) . In Figs. | we show the behavior of T] eff (R) at j8 = 1.0602, 
which is slightly above the estimated value for j8 c , and at jS = 1.1083, which is slightly above the 

(2) 

estimated value for j8 c .A plateau is visible at small distances when L increases and the extension 
of this plateau gets larger with L, consistently with the fact that finite volume effects are becoming 
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Figure 3: rj eff versus R at j8 = 1.0602 (left) on lattices withL = 96, 128, 256,384,512,640,768,1024 and at 
]8 = 1.1083 (right) on lattices with L = 64, 128,256,384,512, 1024. 

less important. The plateau value of T] e ff is about 0.225 at j8=1.0602, i.e. near the first transition, 
and about 0. 16 near the second transition. These values are not far from the expected ones (1/4 and 
4/5 2 , respectively). The determination of T] e ff at j3 c ^ and f}^ is in progress [§]. 

In conclusion, we have determined the critical couplings of the 2D Z(5) vector model and given 
a rough estimate of the critical index v\ near the transitions. Our findings support the standard 
scenario of three phases: disordered, massless or BKT and ordered. In a recent work [ |l0| ] it is 
claimed that the phase transition at J3 C is not a standard BKT phase transition. We will comment 
on this point in Ref. [p]. 
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